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We explore the effects of an external magnetic field in the structure and magnitude of the diquark 
condensate in a three massless quark flavor color superconductor. It is shown that the long-range 
component B of the external magnetic field that penetrates the color-flavor locked (CFL) phase 
modifies its gap structure, producing a new phase of lower symmetry. Our analysis is performed 
within an NJL effective field theory with four-fermion interactions, inspired by one-gluon exchange. 
Using Ritus' method, we compute the quark propagators in the presence of a background magnetic 
field, and derive the gap equations for arbitrary values of the field. An analytical solution is found 
for strong magnetic fields. A main outcome of our study is that the B field tends to strengthen the 
gaps that get contributions from pairs of Q-charged quarks. These gaps are enhanced through the 
field- dependent density of states of the Q-charged quarks on the Fermi surface. Our considerations 
are relevant for the study of highly magnetized compact stars. 

PACS numbers: 12.38.Aw, 12.38.-t, 24.85,+p 

I. INTRODUCTION 

This article presents in full detail the study of the effects of a strong external magnetic field in the gap structure 
of a color superconductor for a system of three massless quark flavors, the main results of which were summarized in 
Rcf. Q. 

Our present knowledge of QCD at high baryonic density indicates that quark matter might be in a color supercon- 
ducting phase (for reviews, see 0)- This interesting possibility has originated a flurry of research. Not only with the 
motivation to deepen our knowledge about QCD, but mainly because of its astrophysical applications. It is very likely 
that quark matter occupies the inner regions of compact stars, or that even quark stars exist On the other 

hand, it is well-known |5j that strong magnetic fields, as large as B ~ 10 12 — 10 14 G, exist in the surface of neutron 
stars, while in magnetars they are in the range B ~ 10 14 — I0 15 G, and perhaps as high as 10 16 G 0. Moreover, the 
virial theorem allows the field magnitude to reach values as large as I0 18 — 10 19 G. If quark stars are self-bound 
rather than gravitational-bound objects, the previous upper limit that has been obtained by comparing the magnetic 
and gravitational energies could go even higher. 

In this paper wc investigate how a strong magnetic field can modify color superconductivity, with the aim of further 
studying its possible astrophysical implications. We will start by considering three massless quark flavors. In this 
case, it is well established that the ground state of high-dense QCD corresponds to the CFL (color-flavor locked) 
phase 0| . In this phase quarks form spin-zero Cooper pairs in the color- antitriplet, flavor- antitriplet representation. 
One can now ask whether this scenario can change when a magnetic field is switched on. Would the external field 
influence the pairing phenomena? As discussed in Ref. Q and shown here in more detail, the magnetic field can 
drastically affect the condensation phenomenon producing a new color superconducting phase that we call magnetic 
color-flavor locking (MCFL) phase. 

To understand how a magnetic field can affect the color superconducting pairing it is important to recall that 
in spin-zero color superconductivity, although the color condensate has non-zero electric charge, there is a linear 
combination of the photon and a gluon that remains massless 0, • A spin-zero color superconductor may thus be 
penetrated by a long-range remnant "rotated" magnetic field B. Although all the superconducting pairs are neutral 
with respect to this long-range field, a subset of them is formed by quarks of opposite rotated charges Q. One would 
expect that this set of condensates will be strengthened by the penetrating field, since their paired quarks, having 
opposite Q-charges and opposite spins, have parallel (rather than antiparallel) magnetic moments. As will be shown 
in the present paper, our intuitive expectation is confirmed by our explicit calculations. The situation here has some 
resemblance to the magnetic catalysis of chiral symmetry breaking jToL ITU , in the sense that the magnetic field 
strengthens the pair formation. Despite this similarity, the way the field influences the pairing mechanism in the two 
cases is quite different as will be shown below. 

Our results signalize the color superconductor as a very peculiar kind of superconductor. Up to now this is the 
only physical system where magnetism and superconductivity are not at odds with each other, but on the contrary, 
the fermion-fermion pairing, which is the hallmark of superconductivity, is reinforced by the magnetic field. 

The plan of the paper is as follows. In Sec. II we discuss the color-flavor and Dirac structure of the gap in the 
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presence of a constant magnetic field. The mean-field theory in the presence of the rotated magnetic field is introduced 
in Sec. III. Special attention is devoted to obtaining the mean-field effective action in momentum space in the presence 
of the external field and to the application of the Ritus transformation ^2] to this case. In Sec. IV we find the gap 
equations and solutions in the presence of the B field. Starting from the general gap equations for arbitrary B values, 
we analytically find the solutions for the limiting cases of zero and strong fields (eB > /i 2 ). We find that at strong 
magnetic fields the gaps that get contributions from pairs of Q-charged quarks are enhanced by the field. In Sec. V 
we sketch the derivation of the low-energy effective field theory of the Goldstone modes in the MCFL phase. The 
main outcomes of the paper and discussion of the results are stated in Sec. VI. 

II. THE MCFL PHASE 

As known, a degenerate quark system becomes unstable under any quark-quark attractive interaction at the Fermi 
surface. This instability leads to color-superconducting Cooper pairing. Depending on how large the density is, there 
are different choices for the interaction. For three different quark flavors, at densities higher than the strange quark 
mass, the quarks can be considered effectively massless. This system of very dense massless quarks can form a favored 
state with spin-zero Cooper pairs in the color-antitriplet, flavor-antitriplet representation (the CFL phase) In 
compact stars where the macroscopic volume of quark matter must be an electrically neutral color singlet, this phase 
is preferred at high enough densities |l3j . 

An important feature of spin-zero color superconductivity is that although the color condensate has non-zero electric 
charge, there is a linear combination of the photon and a gluon G* that remains massless 0,13; 

A lt = cos0A M -aneC^, (1) 

while the orthogonal combination G® — sin 9A fl + cos 9 G* is massive. In the CFL phase the mixing angle 9 is 
sufficiently small (sin# ~ e/g ~ 1/40). Thus, the penetrating field in the color superconductor is mostly formed by 
the photon with only a small gluon admixture. 

The unbroken U(l) group corresponding to the long-range rotated photon (i.e. the U(l) c .m.) is generated, in 
flavor-color space, by Q = Q x 1 — lxQ, where Q is the electromagnetic charge operator. We use the conventions 
Q = — Ag/VJi, where As is the 8th Gell-Mann matrix. Thus our flavor-space ordering is (s, d, u). In the 9-dimensional 
flavor-color representation that we will use in this paper (the color indexes we are using are (l,2,3)=(b,g,r)), the Q 
charges of the different quarks, in units of e = e cos 9, are 
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Although the interaction of an external field with dense quark matter has been investigated by several authors 
0> the effect of the penetrating B field on the superconducting gap structure was not in the scope of these 
studies. However, as shown in our previous paper 1], the B field can change the gap structure and lead to a new 
superconducting phase. To understand how this can occur notice that due to the coupling of the charged quarks 
with the external B field, the color-flavor space is augmented by the charge operator Q, and consequently the order 
parameter of the CFL splits in new independent pieces. In this paper we show how the color superconducting pairing 
is modified in the presence of a magnetic field. 

Let us introduce the rotated-charge projectors 

fi = diag(l,l,0,l,l,0,0,0,l) , (3) 

Q+ = diag(0, 0,0, 0,0,0, 1,1,0) , (4) 

n_ = diag(0,0,l,0,0,l,0,0,0) , (5) 

satisfying 



Cl n £l v i = Sr^tilq, 77,7/ = 0, +, - 



(6) 
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n a + n+ + ^_ = 1 , (7) 

In terms of f2+ and fl- the rotated charge operator takes the form 

Q = 'Af = °+ - Q - ■ ( 8 ) 

77=0, ± 

Using the projectors ©-© we can express the column fermion field in the 9-dimensional color-flavor representation 
(si,S2,S3,di,d2,d 3 ,ui,u 2 ,u 3 ) as the sum 

^ = ip(o) + i>{+) + V'(-) > ( 9 ) 
where subindexes (0)-, (+/— ) indicate fields of zero, positive and negative rotated charge respectively defined by 

-0(0) = , ^(+) = n+i> , ip(~) = &-i> ■ (io) 

The interaction of massless quarks with the external rotated-magnetic field is described by the Lagrangian density 
term 

L e ™ arks =^{XL^)i> , (11) 

with 

n M = id„ + eQA„ . (12) 
Given that the quarks have different rotated charges, the Lagrangian density l|llf> naturally splits in three terms 

L^Zrks = ^ {0 )(id^(o) + ? (+) (id M + el^ {+ ) - eA^(-) ■ (13) 

There is no reason to expect that the color-flavor structure of the condensate will be the same as the CFL one. 
Instead, one would think that the charge-separation of the fields enforced by the external magnetic field should be 
somehow reflected in the color-condensate too. Below we will show that this is indeed the case. 

At this point we need to define a working model where we can investigate the effects of the external magnetic field 
in a system of massless quarks of three different flavors. With this aim, all our calculations will be done in the context 
of a Nambu-Jona-Lasinio (NJL) four-fermion interaction model, abstracted from one-gluon exchange interactions in 
QCD 8J. Although this simple model disregards the effect of the B-ficld on the gluon dynamics and assumes the same 
NJL couplings for the system with and without magnetic field, it keeps the main attributes of the theory, providing 
the correct qualitative physics. 

In order to find a reasonable ansatz for the gap in the presence of a magnetic field, we will start considering the 
system without magnetic field and will find the most general structure allowed for the color superconducting (CS) 
condensate. The corresponding general structure in the presence of the magnetic field will be a particular case of 
the one at zero field, because the two structures have to coincide in the limit of zero field and because the field 
reduces the symmetry of the original theory to a subgroup of the original group. We underline that such a general 
structure is found using only very general symmetry arguments. The gap that minimizes the free energy is just one 
of many possible particular condensate patterns contained within the general structure: precisely the one that retains 
the highest degree of symmetry. This was the guidance principle that led to the CFL ansatz y| in the case without 
magnetic field, and it will be the same principle that will ultimately take us, when a magnetic field is present, to a 
new ansatz: the MCFL one 1]. As seen below, similarly to CFL, in the MCFL phase quarks of all three colors and 
all three flavors pair. The main difference is that the condensate distinguishes between gap parameters that can get 
contributions from rotated charged quarks and those that only get contributions from neutral quarks. 

Let us recall that the gap matrix is formed by diquark condensates (q T Oq), where q T is the transposed field 
operator, and O — Ooirac <S) ® flavor (S) O color is an operator acting through a direct product of the Dirac, flavor and 
color spaces. Fermi statistics constraints O to be totally antisymmetric, i.e. O t = —O. 
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We are interested in a spin-zero condensate and will restrict our attention to the spin-zero channel that is anti- 
symmetric in Dirac indices. Hence, the color-flavor structure of the gap should be symmetric under simultaneous 
exchange of both color and flavor indices. The Dirac structure is then given as C75, where C = 17270 is the matrix 
of charge conjugation. 

The one-gluon exchange interactions are attractive in the channels that are antisymmetric under the exchange of 
color. This fact, combined with the constraint imposed by the Fermi statistics implies that the dominant contribution 
to the gap comes from terms which are antisymmetric under the exchange of color indices and antisymmetric under 
the exchange of flavor indices. This is the color- antitriplet, flavor-antitriplet condensate. It can be shown |l5j that the 
antisymmetric gap can be specified in general by only three gap parameters. Color symmetric channels are repulsive, 
hence subdominant, but once the color- antitriplet, flavor-antitriplet condensate is generated, a color sextet, flavor 
sextet condensate is also induced due to mixing 0] . Following a derivation similar to the one carried out in Ref. , 
one can show that the symmetric condensate can be specified in general by only six gap parameters. 

Therefore, the condensate structure can be written in general as 

A = A Sl Si + A S2 S 2 + A S3 S 3 + A S4 S 4 + A S5 S 5 + A Se S 6 + A^Ai + A Aa A 2 + A A3 A 3 , (14) 
where the six symmetric and three antisymmetric independent color-flavor structures are given by 

Si = 2S al 8 bl 6 a S : > 1 S 2 = 25 a2 8 b2 5 l2 5 32 S 3 = 28 a3 5 b3 8 a 8 j3 S 4 = (8 al 8 b2 + 8 a2 8 bl ){8 ll 8 ]2 + 8 i2 S jl ) 

S 5 = (8 al 8 b3 + S a3 S bl )(S a 8 : > 3 + 8 l3 S jl ) S 6 = (8 a2 S b3 + 8 a3 8 b2 ){8 l2 8 l3 + 8 l3 8 32 ) 



Ai — SabiEiji A 2 — £ a b2£ij2 A3 — £ a b3£y3 (15) 

where a, b denote color indices and i, j denote flavor indices. In the case of the NJL theory, since the quark interaction 
is modelled by a point-like four-fermion term, the nine coefficients A in Eq. Q14f) are independent of the momentum. 

As mentioned above, the energetically favored condensate should be the one retaining the highest degree of symmetry 
possible. In the absence of a magnetic field, that corresponds to the so-called color-flavor locked group SU(3)c+l+r- 
Therefore, the original SU(3)c x SU(3)l x SU(3)r x U(1)b x U(l) e . m . symmetry is broken by the color condensate to 
the diagonal subgroup SU(3)c+l+r- This is the well-known CFL phase, which results from Eq. I|14|l if one demands 
that Agj = ... = As 6 = As and A Al = ... = A^ 3 = A^. The CFL structure parameter can be expressed as 
Acfl = As(U + Uq) + A^(U — Uq) with the color-flavor matrices Uq and U defined by 



Ur 



2 ^ 2 
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Let us consider now the situation with magnetic field. A background magnetic field explicitly breaks the flavor 
symmetry of the original theory to SU(3)c x SU(2)l x SU(2)r x U(1)b x U(l)^ x U{l) e . m ., since u quarks have 

different electric charge than d and s quarks. 17(1)^ ^ is related to the existence of an anomaly-free current given 
by a linear combination of the s, d, and u axial currents The subgroup that ensures preserving the largest 

degree of symmetry in the CS gap is SU{2)c+l+r in this case. Hence, the symmetry breaking induced by the color 
condensate in the presence of a magnetic field is: SU(3)c x SU(2) L x SU(2) R x U(1)b x U^\1)a x U(l) e . m , — > 
SU(2)c+l+r x U(l) e m ■ Notice that this symmetry breaking preserves the rotated electromagnetism U(l) c . 

The SU(2)c+l+r, symmetry requires the invariance of the gap under simultaneous flavor (1 <-* 2) and color (1 <-* 2) 
exchanges. This implies the following equalities among the gap parameters in Eq. I|14f> 

As- = A Sl = A S2 , Af = A Ss = A Se , Af = A Al = A A2 , (18) 
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Based on the above considerations, the gap structure in the presence of a magnetic B field takes the form 
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This is the so called MCFL ansatz. In going from l|14|) to l|19|l we introduced the more convenient notation 
As" = As 3 , As = As 4 , Aa = A^ 3 . At weak magnetic field, B, the CFL and MCFL phases should be continuously 
connected. In the limit of zero field the gap parameters ought to satisfy A B = Aa and Af = As = A$' = As". 

To trace back the physical origin of the different gap parameters appearing in the MCFL structure (|19|l we can 
see that despite the overall Q- neutrality of all the pairs, they can be formed either by a pair of neutral or by a 
pair of rotated charged quarks. Any pair formed by Q-charged quarks feel the background field directly through the 
minimal coupling of the quarks in the pair with B. The gap parameters that get contributions from pairs formed by 
Q-charged quarks must be different from the gap parameters getting contributions only from pairs of neutral quarks. 
The gaps ^-a/s are antisymmetric/symmetric combinations of condensates that get contributions from both charged 
and neutral pairs. On the other hand, the gaps A-a/s,S',S" > are antisymmetric/symmetric combinations of condensates 
that only get contributions from pairs of neutral quarks. The only way the field can affect ^a/s,S',S" is indirectly, 
through the system of highly non-linear coupled gap equations. 

The MCFL order parameter (|19|l can be rewritten in a compact way as 
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(20) 



Out of the four independent, symmetric gaps, all of which are due to subleading interactions, one can single out 
the set of three gaps A s , A s and As that get contributions only from pairs of neutral quarks. There is no reason to 
expect that these three symmetric gaps will differ much in magnitude. Hence, in a first approach to the problem, we 
will solve the gap equations assuming As — A s ~ A s . 



III. EFFECTIVE ACTION IN THE PRESENCE OF A ROTATED MAGNETIC FIELD 

In this Section we introduce the mean-field action in the presence of a constant rotated magnetic field in both 
coordinate and momentum spaces. In our derivations we will follow a parallel approach to the one developed for the 
case of zero field in Refs. [la.ll8|. 

As discussed in the previous section, the rotated magnetic field naturally separates the quark fields according to 
their Q charge. The independent mean-field effective action is 

I B @,il>) = I d 4 xd 4 y{^ [^(Q^^o^ 1 ^' f)^(Q)(f) + ^(Q^^l ^^ 1 ^' v)^(q)o(v) 

Q=0,± 

+ 2^(o)c(*)Aj ) ^ ( o)(y)+?(+)c(*)^ ) V'(-)(»)+?(-)oW^ ) %)(y) + /».cO} , (21) 

T 

where ipc(%) = Ci\> \x) is the charge-conjugate spinor of ip(x). 

In (|21|) symbols in parentheses indicate the correspondence to neutral (0), positively (+) or negatively (— ) Q- 
charged fermions. Propagators for fields and conjugated fields will be denoted with the customary ± supra-indexes. 
Explicit expressions of the inverse propagators are 

[G% )0 ]-Hx,y) = [vfd^ivftS^x-y) , (22) 
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[Gf +)0 }-\x, y) = [vfU^ ± M7 °]<5 4 (z y) , (23) 
[Gf_ )0 ]- X (x, y) = [iynH ± - J/) , (24) 



with 



= ify ± cL4 M , (25) 

and the gap matrices are given by 

A+ 0) = [A s + Aa] + (U - N)Q + [A s - A A ]+UO + [Af + A^]+Nft , (26) 

A+ +) = [Af-A^+)im + , (27) 

A+ ) = [Af -A^]Wun_, (28) 
with the color-flavor matrix N defined as: N = ^[S$ + Sq + Ai + A2] . 

A. Effective action at B 7^ in momentum space 

The computation of the field-dependent quark propagators in momentum space can be managed with the use of 
a method, originally developed for charged fermions by Ritus [12J and later on extended to charged vector fields by 
Elizalde, Ferrer and Incera |19|. In Ritus' approach the diagonalization in momentum space of the Green's functions 
of charged fermions in the presence of a background magnetic field is carried out with the help of the eigenfunction 
matrices E p (x). They are the wave functions of the asymptotic states of charged fermions in a uniform magnetic field 
and play the role in the magnetized medium of the usual plane- wave (Fourier) functions e lpx at zero field. This method 
has also been used in the context of chiral mass generation in a magnetic field fill . l20i | . Using the E p (x) functions, we 
transform the propagators (|23H - il24ll to momentum space adequately. Below we describe the basic properties of the 
transformation. 

The transformation functions Ep(x) for positively (+), and negatively (— ) charged fermion fields are obtained as 
the solutions of the field-dependent eigenvalue equation 

(n( ± ). 7 )4 ±) (a:) = 4 ±) (*)(7-P {±) ) , (29) 

with p^> given by 



and 



with eigenfunctions 



p^ = (p ,0,±^2\eB\k,p 3 ) , (30) 

4 ±) w = E^ ) w a w ' ( 31 ) 

(7 

Eg\x) =K w e-^°+P^+P^) D {e(±)) , (32) 
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where -D n( ±) (f?(±)) are the parabolic cylinder functions with argument Q(±} defined by 



g {±) = ^2\eB\( Xl ± P2 /eB) , (33) 

and index «.(±) given by 

eB 1 

n(±) = n {±) {k,a) = k±—~-a- - , n (±) = 0, f , 2, ... (34) 

2\eB\ z 

k = 0, 1, 2, 3, ... is the Landau level, and cr is the spin projection that can take values ±f . The normalization constant 

W„ (±) is 



N n{±) = (4rr|e[B|)V^n (±) ! ■ (35) 
The spin matrices A(ct) in (|3f |l . not to be confused with the gap coefficients, are spin projectors. They are defined 



as 



A(cr) = diag(<5 CT i,<S CT _i, S al , S a -i), a = ±1 , (36) 

and satisfy the following relations 

A(±) f = A(±), A(+) + A(-) = l ! A(±)A(±) = A(±) , A(±)A( T ) = G, (37) 

7 "A(±) = A(±) 7 ", 7 ± A(±) = A(t)7 X • (38) 
In Eq. ij3~5|) the notation 7" = (7°,7 3 ) and 7 1 - — (7 X ,7 2 ) was used. 

Under the E p (x) functions, positively ("0(+)) and negatively (V'(-)) charged fields transform according to 

= ^(04^ ±} W^(±)(P) , (39) 

= ^j0^(±)(P)Ef\ x ) > (4°) 

where E?\x) = 7o (4 ±) (*)) t 7o and £ ^ = £~ / 
One can show that 

[7„(n c+)M ± ^ )]4 +) (x) = 4 +) (z)M?i +) ± M^o)] , (41) 

and 

[ 7M (n ( _ )M ± /x^o)]^^^) = El-\x)b M ft-) ± M^o)] • (42) 

Since the charge conjugate of a positively (negatively) charged field is a negatively (positively) charged field, the 
charge conjugate fields transform as 

^(+)o(a) = ^-^l4~H^(+)c(p), (43) 
^-)c^)=^^K +) (^(-)ciP) ■ (44) 
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This transformation dictates what fields should form the components of positively and negatively charged Nambu- 
Gorkov fields. 

In terms of Nambu-Gorkov fields ^, the effective action in the presence of magnetic field B takes the form 

I B {1>,4>) = J d A xd i y^{x)S- 1 {x,y)^{y) 



(45) 



where 



«S ( o)(p) = 



(46) 
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The Nambu-Gorkov fermion fields corresponding to the different Q charges are defined by 



*o = ( ^ ) 



for the neutral field, 



^(+) 
V>(-)c 



for the positive field, and 



(v>(+)c) 



(48) 



(49) 



(50) 



(51) 



(52) 



for the negative field. 

As noticed before, the positive (negative) Nambu-Gorkov field is formed by the positive (negative) fermion field 
and the charge conjugate of the negative (positive) field. This means that, as it should be, the rotated charge of the 
up and down components of a given Nambu-Gorkov field arc the same. This in turn determines the nature of the 
Q-charges of the fields entering in a given pair, which we know should have zero overall Q charge. 

The bare inverse propagator of the neutral field is 



[ G (o)o] X (P) = 7 iU (?V ± M<W 



(53) 



Here the momentum is the usual p — {po,Pi,P2,Pz) of the case with no background field. 
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For the positively and negatively charged fields the bare inverse propagators are 

[Gf +)0 ]- 1 (p)=^(p^±^ ) , (54) 

[Gf_ )0 ]- 1 (p)=7p(^ ) ±M^o) • (55) 

respectively. 

B. Nambu-Gorkov propagators in a color-flavor rotated basis 

In order to find the propagators in Nambu-Gorkov space one has to invert the matrices l|46 [l -l]48 [l . whose color-flavor 
structure is rather complicated. At this point we find convenient to perform a rotation that can partially diagonalize 
the components of the Nambu-Gorkov propagators in color-flavor space. 

Let us define a new basis for the fermion fields 



where the indices a and i refer to color and flavor respectively, and X A , for A = 1, . . . , 8 are the Gell-Mann matrices, 

while A 9 = yfl. Given that in the CFL case the Nambu-Gorkov quark propagators become color-flavor diagonal 

pl| in the basis we call l|56[) the CFL basis. 

In the CFL basis the quark kinetic term (flip becomes 

Quarks = # V (i9„S AB - e Q^ fl i M ) ^ B , (57) 

where the charge matrix is 

Q AB = -±Tr([\ A ,\ B ]Q) . (58) 

In our conventions, Q = — -^|, and thus Q AB = ~^f Am is expressed in terms of the antisymmetric structure constants 
of SU(3). In the basis l|56|) one can as well find the charge-projector operators 

n^=diag(l,l,l,0,0,0,0 ) l,l), (59) 

and 

n AB = ^[(S A4 S m + S A5 S m + S A6 S m + 6 A7 5 BJ ) ± i(8 AA 5 m - 6 A5 S m + 8 M 5 B1 - S A7 S m )] , (60) 

which obey the conditions Furthermore, one has (fl±) T = so that Q±ipc — (V't)c- 

From Eq. (|56(l one can deduce the transformation law of the gap matrix, which now reads 

a as = \(\ A )lA l LK ■ (61) 

After taking into account the approximation As ~ Ag/ ~ As" in the gap matrix l|19|) one obtains the gap parameters 
in the new basis to be 

(62a) 
(62b) 

(62c) 

f) , (62d) 
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while all the remaining elements of A AB are zero. 

In the CFL basis the Nambu-Gorkov effective action (|45|l can be written as 

«$,1>) = J ^*£(P)P ( 0)(P)]^*?(P)^ (63) 

As follows we find convenient to express all propagators in terms of energy projectors. For completion, in Appendix 
E|we present the details on how they are generalized for charged particles in the presence of an external magnetic 
field. 

For the charged fields the Nambu-Gorkov matrix reads 



0,,=," (%...»;;;> %>..»;;;;,, m 



where the indexes A,B take values 4,5,6,7 only and the subindices 11,12,21,22 specify the entry of the Nambu- 
Gorkov matrix. Notice that S^(p) depends on the Minkowski momenta p^ = (po,P^) defined in Eq. (|3(JI) . To 
avoid confusion, we will always denote with caligraphic letters the Nambu-Gorkov matrixes, and with capital letters 
the elements of the matrixes. 

For massless quarks, the elements of the Nambu-Gorkov matrix for the positive charged quarks are 

„ H+) = Af +) (P)7°(P0TM±|P (+) I) Af +) (P)7°(P0TMT|P (+) |) 

(+)ll/22lP j P ^(M-|pW|) 2 -(A5-Af)2 + p2_ (M+ |p( + )| )2 „ (A B_ Af)2 > ^) 

q H+h _ f Ar +) (P)(A|-Ag) A+ +) (p)(Af-AjQ \ 

i+hAP ] ~ 75 l^-(^-|pW|) 2 -(Af-Af)2 + ^- (A1 + |p( + )| )2 -(Af-Af)^ ' [bbb) 



q m+K - J A+ +) ft)(Af-Agr A r+) (p)(Ag-Ag) ■ 

i+h2[P j " 75< i^-(/,-|p(+)|) 2 -(Af-Af)2 + ^_ (M + |p( + )| )2 _ (A B_ Af)2 ^ ^ 



Following the definition of the momenta (J20J), we have | p^- 1 1 = y2|ei?|fc+p§. We use the symbol A^ (p) to denote 

positive/negative energy projectors of positively charged quasiquarks in the presence of the external field, as defined 
in Appendix lAl 

The propagator for the negatively charged fields keeps the same structure, with the only change p^ — > p^~\ Its 
off-diagonal term satisfies the equation, 

S(-) 21 (x, x) = S {+hl (x, x) , (66) 

since 

y-k d 4 p ± j A-(±)(P) (Af - Af) 

S(± hl (X, X) = ^J^E p (x) 7 5 | p2_ (A1 _|p( ± )| )2 _ (A B_ A B )2 



A+ (p) (Af - A 



^-(^+|p (±) l) 2 -(Af -Af)2 
so after integrating in pi and using the relation 



E±{x) , (67) 



eB 

dp2D n {pt±))D n >(pr±- ) ) = nW2-K—== S nn > (68) 

V2eB 



we obtain 



eB^ f dp dp 3 B B ( 1 
S(±hAX ' X) = " A ^ 75 {p§-(M-|p|) 2 -(A^-Af) 2 



p2-( M+ |p|)2_ (A B_ Af)2 J ' 



(69) 
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where |p| stands for the common modulus of the three-momentum of positively and negatively charged quasiparticles, 
since |p| = |p + | = |p~|. From now on we will use the short notation |p| for the modulus the three momentum of 
charged quasiparticles. 

The propagator for the neutral fields in the CFL basis with indexes A, B = 1, 2, ...9 is given by 



cAB 



(p) 



\ 



o 










s, 






33 



0000 
0000 



#)(P) 00 



0000 
0000 
0000 
0000 











89 (p) 



00 00«S ( 8 8 } (p) s\ 
00 S™(p)S$(p) J 



where each component is a 2 x 2 Nambu-Gorkov matrix 



?AB 



One finds that 



o22 



<?33 



<Sfm (p) 



(0) 2 1 (?) S (0) 2 

with Nambu-Gorkov components given by 

A ± (p) 7 °(poTM±IpI) , AT(p) 7 0( PoTMT |p|) 



Pi 



oil 



(0)2 



l (P) = 75 



A"(p) (A s - 



A A ) 



pl - (^) 2 

A+(p) (A S -A A ) 



P 2 , 



"75 



(e£) 2 

f A+(p) (A s -A A y 



P 2 , 



+ 



A"(p) (A s -A^) 4 



P 2 , 



(^) 2 



with 



(70) 



(71) 



(72a) 
(72b) 
(72c) 



e ° = V(m-|pI) 2 + (Aa-a s )2 , f p = ^AJ+W 

The remaining block matrix for the CFL coefficients A = 8, 9 in (|70|l 



(A A -A S ) 2 



(73) 



' c88 c89 
°(0) °(0) 
<?98 <?99 

v (0)°(0). 



(74) 



has matrix elements which are 2x2 matrixes in the Nambu-Gorkov space. Given that the symmetric gap As comes 
from subdominant channels, we can assume, as in the CFL case, that Ag <C A^. Using this approximation the 
Nambu-Gorkov components of each of the matrix element of l|74|l can be expressed as 



S, 



AB 
(0)u 



qAB 
b (0) 



A ± ( P h°(poTfi±\p\)A AB , AT(p) 7 "(p oT/iT | p |)i AB 




Pl (z a P ) 2 ) (Pl- (4) 2 ) 
A"(p) B AS 



(p^-(^) 2 )(p^-(^) 2 ) 
A+(p) B AB 



(Pl (<# 2 ) (P 2 - (e£) 2 ) (p 2 - (e«) 2 ) (p 2 - (^ 
A+(p)(B AB )* A-(p)(B AB )* 



{Pl ~ (e£) 2 ) (p 2 , - (e£) 2 ) [pl ~ (e a P ) 2 ) (p| - (e£) 5 



(75a) 
(75b) 
(75c) 



where 



^" = J(n- |p|) 2 + (A a/fc ) 2 , = 



v /( / ,+ |p|)2 + (A a/6 ) 2 



(76) 
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with 



1 



A 2 /h = - (4A% +A\± A A ^A\ + 8A 2 , j . (77) 
and we have defined for convenience 

A N = A* + Af . (78) 

The CFL submatrices A and B read 

/r 2 -|(A A -A w ) 2 -|(A A + 2A,v) 2 f A A (A A - A N ) \ 

I ^A A (A A -A N ) r 2 -l(A A -A N f-^(A A ~4A N f ' 



B = (A A - 4A N ) + IA% (A A + 2A N ) & (A A - A N ) (r 2 - 2A 2 N ) \ 
^ ^ (A A - A N ) (r 2 - 2A 2 N ) 2 £(A A + 2A N ) + ^(A A ~4A N )J ' 

where we used the following shorthand r 2 = p\ — (|p| — n) 2 . The matrices A and B are the corresponding antigap 
quantities. They are obtained simply by replacing r 2 — * f 2 = p^ — (|p| + ji) . 



IV. GAP EQUATIONS AND SOLUTIONS IN THE PRESENCE OF A ROTATED MAGNETIC FIELD 

The evaluation of the different fermionic gaps can be done by writing the Schwinger-Dyson equations. In coordinate 
space the QCD gap equation reads 

A + (x,y)=i^Xl^S 2 i(x,y)rX B D^(x,y) , (81) 

where, for simplicity, we have omitted explicit color and flavor indices in the gap and fermion propagator. Here 
is the gluon propagator. 

We leave for a future project the study of Eq. I|81|l in the high density /weak coupling limit of QCD in the presence 
of a magnetic field. It should be noted, though, that gluons can also feel the presence of a rotated magnetic field 
because some of the gluons carry Q-charges, and thus they minimally couple to the external £?-field. Debye and 
Meissner masses, as well as Landau damping, should then be affected by B. Hence, before studying the QCD gap 
equation l|81|l in the presence of a magnetic field, it is necessary to study the behavior of the gluon propagator for 
different energy scales, so as to properly assess its effect in the gap equation. 

In the NJL model abstracted from one-gluon exchange that we are using in this paper the gap equation can be 
obtained from Eq. I|81l) simply by substituting the gluon propagator by 

tf? (x, y) = i 5 AB 6^ (x-y). (82) 

and using the quark propagator in the presence of the rotated magnetic field found in Section IlIII 

The NJL model is characterized by a coupling constant g and an ultraviolet cutoff A. The ultraviolet cutoff should 
be much larger than any of the typical energy scales of the system, such as the chemical potential fj, and the magnetic 

energies \J £B. In other studies of color superconductivity within the NJL model the values of g and A were chosen so 
to match some QCD vacuum properties, and consequently hoping to get also the correct approximated quantitative 
results for the gaps. We will follow the same philosophy here, however, noticing again that this completely ignores 
the effect of the magnetic field on the gluon dynamics. 



A. Gap equations for arbitrary values of the magnetic field 

Using the transformation rules given in Eq. l]5(:iH one can re-write the NJL gap equation in the CFL basis 

(A+f c = £ {Tr (\ S \C\ c \ b ) i T r (A^) Tr (X c X B ) \ j^Sf^ (*, x)^ . (83) 

B,C=1 
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After computing the above 17(3) traces, taking into account the particular CFL-basis structure of the propagator, 
where = S%f = S^f, and keeping in mind that the relation 1|66|) is held, we can write down more explicitly the 
gap equations for every element of A AB . They read 

(A+) 11 = -z^r {-^(oV (x, x) + \ (Sf$ )2i (x, x) + 2V2Sf* )2i (x, x) + 2S% 21 (x, xj) } 7/l , (84a) 

(A + ) 44 = -^l" {-~ (2S?$ hl (x,x) + V2S^Jx,x) - 2S™Jx,x)) - |s (+)21 (x, x)} 7m , (84b) 

(A+) 88 = -i^V {s$ )ai (x,x) + I (Sf* )21 (x,x) - 2V2Sf° )21 (x,x) + 2S^Jx,x)) - |s (+)al (x, x) } 7 „ . (84c) 

(A+) " = ^IF 7 " { 25 (<U^) + j^ifo*) + . (84d) 

( A+ ) 89 = -i^y {sS^Cs.ar)- (^(s.x) + VS^^a:)) - 2S (+)ai (x >iB )} 7 ^ . (84e) 
From Eqs. lEU one can immediately deduce 

A A = ~(5A n - A 88 - A 99 ) , (85a) 

A s = - (A n + A 88 + A 99 ) , (85b) 
6 

A N = Af + Af = i (A n - 5A 88 + 4A 99 ) , (85c) 



Ac = Af - A A = A 44 ■ (85d) 



Therefore, 



y S^ )91 (x, x)+ 1 - (S$ hl (x, x) + 6V2Sf* )21 (x, x) + 4S^ )21 (x, x)) } 7m , (86b) 
(Af )+ = -i^-f {S#)jx,x) ~ 5Sf* hi (x,x)+±Sf* hi (x,x) + l2S i+)21 (x,x)}^ , (86c) 

2 

«) + = "^7 P {S$ )n {x,x) + 7S$ )al (x,x) + 6V25 8 9 )2l (x,ar) - 8S^ 2i (x,x) + 2AS (+hl {x,x)) 7m , (86d) 

Using the explicit expressions of the propagators given in Sec. IIII 51 we can write the gap equations in momentum 
space. In the limit Ag <§C A^, and dropping the contribution of antiparticles, Eqs. (|86a|l . I|86b(l . I|86c|l . and (|86d|> are 
respectively given by 

_ . ff 2 f rf 4 g f A A C a + M Cb — M \ 

S 9A*J (27r)H«o 2 -(|q|-M) 2 -Ai 9o 2 - (|q| - M ) 2 - A* gg - (|q| - - Af / ' ^ 
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_ . g 2 f d A q f 17 A A 1 7C a - 2M 1 7C fe + 2M ] 

A *9A 2 7 (2^\2g2_ ( | q |_ M)2 _ A 2 4 + 2 (? 2_ ( | q |_ / ,)2_A2 + 2(? 2_ ( | q |_ M) 2_ A 2j . l«< J 

A B = _•_£__ / / A f F a +9£>„ F b + 9£> b 

S ? 36A 2 [J (27T)4 \ gg _ (| q | _ M) 2 _ + gg _ ( | q | _ M) 2 _ A 2 + gg _ (| q | _ ^2 A 2 

+ UeB y (2,)3 g0 2 -(i q i-,) 2 -A 2 c ) • ( 87c ) 

A *36A 2 V (2^ g 2 - (|q| - - A\ ql (|q| - M ) 2 - A 2 + gg - (|q| - M ) 2 - A 2 

where A 2 ^ h was defined in Eq. (|77jl . and 

_ A A / 4A 2 , + Aj \ 2Af v 

Ca/fc -~rWATTMfJ ' m -7attmT (87e) 

Fa/fc = A^6A, ± A 2 ^ OA v A., > (8?f) 



2^A 2 4 + 8A 



O a/t = A w (l± /a ., Aa „ ,„• 1 ■ 



V&a + 8A 



These are highly coupled non-linear equations, that can only be studied numerically. However, one can find a 
limiting situation where an analytical solution can be found. It corresponds to the situation of moderately high B 
fields (fi 2 <eB< A 2 ). 

B. Zero magnetic field limit 

For completion we show here that in the limit of zero magnetic field the gap equations l|87a|) - (|87d|) reduce, as it 
should be, to the CFL gap equations. 

In the B — > limit, one should have A^, g — > &A/S- Using 

A^ = Ayi » Af = A s , (88) 
in the definitions (|77|l and (|87e|) -( |57g| ), they take a more simplified form given by 

A 2 «4Ai, A 2 ^, C a *^A A , C b n-±A A , 
F a « -4A.4, F b » -A A) D„ « ^A/i, D b « ^A A , M » ^A A . (89) 

Making use of (|88|l and l|89l) in the right hand side of Eqs. (|87al87d() and taking into account that when B — ► the 
Landau level summation should be replaced by an integral following the prescription 



eB£/(2eB»)- / ^/(d.) 00) 



2tt 
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we can easily see that both Eqs. (|87a() and l|87c|) reduce to the CFL gap equation for the symmetric gap 



f A \ 2A A 2A A 

*9A 2 J (2^)4 \ ql (|q| - M )2 _ A 2 A q 2 - (|q| - 4A 2 



(91a) 



while both Eqs. (|87b|) and (|87d|) reduce to the CFL gap equation for the antisymmetric gap 



A, = 



^ 2 r d 4 q 



8A, 



4A, 



9A 2 ./ (27r)Ug 2 -(lq|-M) 2 -Ai 9 2 " (IqI " M) 2 - 4A 2 



(91b) 



Therefore, starting from the Eqs. H87a[l - I)87djl we recover the correct CFL gap equations in the NJL theory. The 
quark propagators in the CFL case can be straightforwardly found, (see for example Ref. H2), assuming that 
As <C Aa- The solution reads (S = A — [m) 



A^ fl ~ exp 



and 



3A 2 7r 2 
2<?V 



n 2 ii 2 

aCFL xCFL 9 M , „ 

As Aa 9A^ M2 



(92) 



(93) 



C. Strong magnetic field limit 

Let us now consider the MCFL gap equations for the case eB > n 2 . Taking into account that the main contribution 
to the gap comes from quark energies near the Fermi level, it follows that in the strong field region the leading 
contribution comes from the lowest Landau level (LLL). 

Assuming A^ > Af , A^, and A^ > A s , the definitions (J77JI, (|S7e|) -( |87g| l reduce to 



Fn 



■(3- 



A 2 



2(A 
F b 



B\2 
A) 



"(3 



Al « 2(A*) 2 , 



V2 B 



2 / — ^i 3 ~ u v" 2 7 — ^ 1 — w — ^ 7 — w — ^ 7 2 
In this approximation the gap equations decouple and the equation l|87d|) for A^ becomes 



a 



a/b 



Aa 
2 



A B 



M: 



(94) 



g 

3A 2 



d 3 g 



(2tt)* 



(| q |- M ) 2 + 2(A^) 2 



ff 2 eB 



^A 



3A2 ( 2 -) 2 v /(| (73 |- M ) 2 + (Af) 2 



(95) 



Its solution reads 



A^ ~ 2a/J/I exp ( 



5 2 U 2 



(96) 



The exponent in (|96|l has the typical BCS form exp [ 1/NG ], where G = g 2 /3A 2 is the characteristic effective 
coupling constant of the 3 channel [23] , and N represents the total density of states at the Fermi surface of the four 
quasiquarks of a single chirality contributing to the antisymmetric gap parameter A^ . If no magnetic field is present 
(CFL case), N = 4A M , where = fi 2 /2n 2 is the density of states of one quasiquark. At nonzero magnetic field 
the density splits in two terms N = 2N fi + 2N^, because two of the four quasiquarks are neutral, thus they have 
density Nu, while two are charged and their density of states depends on the magnetic field. In the zero Landau 
level the density of states of a charged quasiquark is Ng = eB/2ir 2 . Notice that the gap parameter l|96|l increases 
with the magnetic field. Comparing Eas. (|92|) and (|96|l one can see that A^ > A^ FL if eB > /i 2 , as assumed in the 
derivation of Ea. (|95|l . This means that the field not just changes the gap structure, but it also enhances the gaps that 
get contributions from pairs formed by charged quarks. 
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As mentioned in the Introduction, although this situation has some similarity with the magnetic catalysis of chiral 
symmetry breaking |l0j ; the way the field influences the pairing mechanism in the two cases is quite different. The 
particles participating in the chiral condensate are near the surface of the Dirac sea. The effect of a magnetic field 
there is to effectively reduce the dimension of the particles at the lowest Landau level, which in turn strengthens their 
effective coupling, catalyzing the chiral condensate. Color superconductivity, on the other hand, involves quarks near 
the Fermi surface, with a pairing dynamics that is already (1 + l)-dimensional. Therefore, the B field does not yield 
further dimensional reduction of the pairing dynamics near the Fermi surface and hence the lowest Landau level does 
not have a special significance here. Nevertheless, the field modifies the density of states of the Q-charged quarks, and 
it is through this effect, as shown in Eq. Q96[l . that the pairing of the charged particles is reinforced by a penetrating 
strong magnetic field. 

In the strong field approximation the remaining gap equations (|87af) - (|87cf) are respectively given by 

g 2 f d 3 q / A A A A \ 

S ^i A (2^lV(|q|- W + *T ,/(|q| - M ) 2 + 2 (Af P ' ^ 



A 



and 



" - 4A 2 A (2.)3 V 9 V(|q| — n) 2 + A 2 A 9 y^p^T ! ' 

9 2 f d 3 q A B A g 2 eB f A dq 3 A A 



^ " " 6A2 A C*0 8 yj (|q| - ^) 2 + 2(Af ) 2 + L ^ yJ^-.fH^) 2 ' 

The corresponding solutions are 

A s ^A(i_ ^ (100) 
17 l + y 

, 21 

A^ v ^exp h - + T7 - (TT - y ], (101) 

Af * (j^ ~ \)A% (102) 

where x = g 2 n 2 / A 2 tt 2 , and ?/ = eB/fi 2 . As expected, the solutions (|96|l . (|l()()(l - l|102(l are consistent with the initial 
assumptions A A 3> Af , A^, and A^ 3> As if the field satisfies eB > /i 2 . When the field increases within this strong 
field region A A grows and A^ decreases, and they clearly split. Even for fields just slightly larger than /i 2 the results 
are reliable. For example, for y = 3/2 and x ~ 0.3 [24|, one finds A^ ~ 0.2 Af for y = 3/2, while for x ~ 1 then 
A A ~ 0.5 Af . For fi ~ 350 - 400 MeV, one can estimate the field strength for y = 3/2 to be eB ~ (1.2 - 1.6) • 10 18 G. 

V. EFFECTIVE FIELD THEORY FOR THE GOLDSTONE BOSONS OF THE MCFL PHASE 

In the absence of a magnetic field, three-flavor massless quark matter at high baryonic density is in the CFL phase on 
which diquark condensates lock the SU(3) color and SU(3) flavor transformations thereby breaking both symmetries 
according to the following symmetry breaking pattern 

SU(3) C x SU(3) L x 517(3)* x U(l) B -> SU(3) c +l+r ■ (103) 

Nine Goldstone bosons appear as the result of the Anderson- Iliggs mechanism. One is a singlet, scalar mode, associated 
to the breaking of the baryonic symmetry, and the remaining octet is associated to the axial SU {3) A group, just like 
the octet of mesons in vacuum. 

Strictly speaking, once electromagnetic effects are taken into account, the flavor symmetries SU(3)l x SU(3)r of 
the original three-flavor massless quark matter are reduced to SU(2)l x SU(2)r x U{l) { A ] with U(l) { £ c SU(3) A , 
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since the electromagnetic charge of d and s quarks is different from that of the u quark. However, because the 
electromagnetic structure constant a e .m. is so small, this effect is really tiny, a small perturbation, and one can still 
consider the original flavor as an approximated symmetry of the theory. 

On the other hand, we know that in the CFL phase a so-called "rotated" electromagnetism remains, allowing the 
penetration of the color superconductor by a conventional magnetic field in the form of a "rotated" magnetic field. If 
the penetrating magnetic field is strong enough, the effect of electromagnetism cannot be treated anymore as a small 
perturbation since, as shown in previous sections, it affects the pairing phenomena and leads to the development of 
the MCFL phase, whose symmetry breaking pattern is 

SU(3) C x SU(2) L x SU(2) R x £7(1)^ x U(l) B x [7(l) e . m . - SU(2) C+L+R x U(l) c m . (104) 

As mentioned in Section ^ the locked SU (2) group corresponds to the maximal unbroken symmetry, such that it 
maximizes the condensation energy. The counting of broken generators here, after taking into account the Anderson- 
Higgs mechanism, tells us that there are only five Goldstone bosons, all of which are neutral with respect to the 
rotated charge. As in the CFL case, one is associated to the breaking of the baryon symmetry; three Goldstone 
bosons are associated to the breaking of SU(2)a, and another one is associated to the breaking of 11(1)^ . Therefore, 
apart from modifying the structure and magnitude of the gap, the applied strong magnetic field does also affect the 
low energy properties of the color superconductor, since it reduces the number of Goldstone bosons from nine to five, 
none of which is charged. 

Following an approach similar to the one discussed in [2lJ and [25| we can derive the low-energy effective theory 
for the Goldstone bosons in the MCFL phase. Here we will just sketch the main steps leaving detailed calculations 
for a future work. The main idea is to introduce left and right coset fields 

X ta ~ e ijk e abc (t/j^tp c L k y , Y ia ~ e tlk e abc {i; b liji R k y (105) 

where a, 6, c denote flavor indices, i,j, k denote color indices, and L/R denote left/right chirality, respectively; which 
transform under the global transformations of the original group as 

X -> U L XUl , Y -> U R YUl . (106) 

where Uc is a color transformation, while U L G SU{2) L x U(l) { ^ and U R € SU(2) R x U(l)<£> are left and right flavor 
transformations respectively. We are ignoring in our discussion the baryon symmetry, but its Goldstone mode can be 
incorporated in a similar way. 

Then, keeping in mind that the Goldstone modes are related to motions along the degenerate minima of the 
effective potential, one can obtain the effective action for the Goldstone modes just by factoring out the modulus of 
the condensate and allowing the phases of X and Y to depend on time and space. 

At this point it is convenient to introduce the flavor singlet 

E = Xyt = exp ( <-?- + # J , $ = ^cr A , A =1,2, 3. (107) 

V U,B J 

where a a are the SU(2) generators and the fields 4>a are associated to the breaking of SU(2)a, while <fio is associated 
to the breaking of U(l)^ . 

Taking into account that the effective low energy Lagrangian of the Goldstone bosons in the MCFL should be 
symmetrical under the rotations of the original group SU(3)c x SU(2)l x SU(2) r x 17(1)^ , we can write it down 
up to leading order in the derivatives as 

C = tlf. (Tr (doXdo^) + (vigf + «f 5 jf ) Tr (fl^tf)) . (108) 

Notice that we introduced different longitudinal and transverse velocities. This should be done because besides 
the usual Lorentz symmetry breaking proper of the theory at finite density, the strong magnetic field induces and 
extra symmetry reduction, since only the spacial SO(2) rotations in the plane perpendicular to the magnetic field 
are allowed. The constants f 7r _ R , v± and v\\ should be computed from the microscopic theory. We leave such a 
computation for a future project. 

It is remarkable that although the external magnetic field does not explicitly enter in the low-energy Lagrangian 
(|108fl . since it does not couple minimally to the MCFL Goldstone bosons because these are all neutral with respect 
to the rotated charge, its presence is manifested anyway through the anisotropy of the velocities. 
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The effective field theory for these Goldstone fields is totally analogous to that found for QCD in a strong magnetic 
field in Ref. 01- Were it not for the presence of the extra Goldstone boson associated to the baryon symmetry 
breaking in the MCFL phase, the two low-energy theories, at low and high densities, in an external magnetic field 
would be completely equivalent, and one would expect them to be connected by a crossover. We might think that 
because the existence of the extra Goldstone mode, this is not possible. However, the situation is similar to the case 
without magnetic field. There also the high and low density effective theories of the Goldstone modes are identical 
with the exception of the extra Goldstone mode in the color superconductor due to baryon symmetry breaking, and 
yet it has been argued 26] that the extra mode in the CFL phase may have a counterpart in the chiral phase, the 
dibaryon state H found in Ref. 27]. Therefore, the question about the possible equivalence between high and low 
density low energy theories in a magnetic field still remains open. 

VI. CONCLUDING REMARKS 

In this paper we investigated the effect of a magnetic field on the color superconducting gap of a dense, three 
massless quark flavor system. Our main result is that the magnetic field leads to the formation of a new color-flavor 
locking phase, that we have called the MCFL phase. We found that the magnetic field affects essentially the gap by 
modifying the density of states of the charged quarks on the Fermi surface. 

In the analytic calculations performed in this paper we assumed a strong field approximation (eB > /i 2 ), which 
corresponds to fields ~ 10 18 G. The only reason we considered such strong fields was to simplify our task so we could 
avoid to sum in Landau levels, and hence could perform analytical calculations. However, we underline that there is 
nothing special about the LLL here, because the main effect of the field is to modify the density of states and that is 
not restricted only to the LLL. On the other hand, there are indications, based on the study of the low energy degrees 
of freedom of the CFL effective field theory in the presence of an external magnetic field [23, that for eB oxt ~ 5 • 10 16 
G the symmetry pattern of the theory qualitatively separates from that of the CFL phase and becomes characteristic 
of the new MCFL phase. These estimates suggest that a distinguishable separation between the gaps and may 
take place already at those magnetic field orders. Of course, a definitive determination of the field strength required 
to separate the two phases can only be obtained through a numerical calculation on which the effect of all the higher 
Landau levels in the gap equations are included. 

We have shown that the MCFL phase has a smaller vector symmetry than the CFL phase and consequently the 
number of Goldstone bosons reduces from 9 to 5. Our zero-temperature results imply that a propagating rotated 
photon with energy less than the lightest charged quark mode cannot scatter, since all the Q-charged quarks acquire 
a gap and all the Nambu-Goldstone bosons are neutral. The anisotropy present in the background of an external 
magnetic field and the existence of charged Goldstone bosons in CFL but not in MCFL indicates a rather different 
low energy physics in these two phases and it should be reflected in the transport properties. In particular, the MCFL 
superconductor is transparent and behaves at T = as an anisotropic dielectric, as opposed to the isotropic dielectric 
behavior of the CFL phase [52, l29|. However, similar to the CFL, the medium will become optically opaque as soon 
as leptons are thermally excited j30fl . 

Although the Q neutrality of the MCFL is guaranteed without having to introduce any electron density due to the 
absence of gapless modes and the Q-neutrality of the diquark condensates, imposing the 3- and 8- color neutralities 
may require the introduction of small "chemical potentials" ^3 and /zs due to the difference between the main gaps 
A A and Af . 

If we introduce quark masses the scenario will be much more complicated though. First, we would have to revise 
the neutrality conditions and see how they are affected. In addition, one would expect that the effect of the magnetic 
field increasing the gaps formed by Q-charged quarks might have some impact in counteracting the stress produced 
by the mass of the strange quark on the corresponding Fermi sphere and perhaps also in the onset of gapless modes, 
along with the problem of the instability of the gluon modes due to imaginary Meissner masses. Considering the 
consequences of incorporating quark masses, together with a careful study of the effects of the magnetic field in the 
low energy physics, in transport properties or in neutrino dynamics will be the subject of future investigations. 

Even if only at a speculative level, it is worth to mention some possible astrophysical consequences of the main 
results of this paper. The equation of state of the color superconductor will be affected by the magnetic field, as 
different quark gaps will be affected in a different way by the field. We do not expect this to be a pronounced effect, 
though, but it might be interesting to see whether it affects the mass-to-radius ratio of the star. Because the low 
energy physics of the CFL and MCFL are so different, depending on the strength of the star's magnetic field there 
may be different signatures in the cooling process of the star. Transport properties, such as viscosities or thermal 
conductivities, will be also affected by the presence of the magnetic field. Finally, the dynamics of the magnetic field 
itself might be very peculiar in quark matter, differing from that in a neutron star, which is commonly believed to be 
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an electromagnetic superconductor. 
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APPENDIX A: PROJECTORS IN THE PRESENCE OF A BACKGROUND MAGNETIC FIELD 

In the absence of a background magnetic field the condensate is usually given in terms of different projectors |16| . 
The chiral and helicity projection operators are given respectively by 

Pn, L - (Al) 



and 



H ± (p) = ^[1± 757CTY • (A2) 



The massless quark free-energy projector for positive and negative energies is given by 

A*(p) = 1±7 f ,? , (A3) 

For massless quarks, the chiral and energy projectors are enough to specify the quark propagators. Furthermore, 
in a NJL theory, where there are only contact interactions, the gap is a momentum independent constant, and the 
Dirac structure of the condensate is particularly simple in the spin zero case, where it simply reduces to C75, where 
C = 27270 is the matrix of charge conjugation. 

In the presence of an external B field the helicity-projection operators (|A2|I and the free-energy projectors (|A3|I are 
not conserved quantities, since they do not commute with the field-dependent Hamiltonian. But their generalization 
in terms of the field dependent momentum operator (|12(1 are. In covariant form, these conserved operators can be 
expressed respectively by 

g ±(n)= 1± ^ g --"^, (A4) 



X ±( n) = 1± ^ n -- u - n ^, (A5) 

where it M is the four-velocity of the center of mass of the many particle system. In the rest frame, the helicity and 
energy projectors can be expressed in momentum space respectively as 

1 -x- - (±) 

5-(±),_s 1 ± 75707 ■ P { \ c\ 

H ± '(p) = = , (A6) 



1 _i_ 

T± /-\ l±7o7'P / A7 x 
A (±)(P) = 2 ■ ^ > 

Similarly to the free case, in the presence of a "rotated" magnetic field the chiral IjAlfl . helicity i|A2() and energy 
(|A7I) projectors commute. 
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